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Summary 

Inconsistent citations in the literature and questions about convergence prompt reexamination of 
Hencky’s classic solution for the large deflections of a clamped, circular isotropic membrane under uni- 
form pressure. This classic solution is observed actually to be for uniform lateral loading because the 
radial component of the pressure acting on the deformed membrane is neglected. An algebraic error in 
Hencky’s solution is corrected, additional terms are retained in the power series to assess convergence, 
and results are obtained for two additional values of Poisson’s ratio. 

To evaluate the importance of the neglected radial component of the applied pressure, the problem 
is reformulated with this component included and is solved, with escalating algebraic complexity, by a 
similar power-series approach. The two solutions agree quite closely for lightly loaded membranes and 
diverge slowly as the load intensifies. Differences in maximum stresses and deflections are substantial 
only when stresses are very high. The more nearly spherical deflection shape of the membrane under 
true pressure loading suggests that a near-parabolic membrane reflector designed on the basis of the 
more convenient Hencky theory would not perform as well as expected. 

In addition, both theories are found to yield closed-form, nonuniform membrane-thickness distribu- 
tions that produce parabolic middle-surface deflections under loading. Both distributions require that 
the circular boundary expand radially in amounts that depend on material and loading parameters. 

Introduction 

Concepts for orbiting inflatable reflectors are of interest primarily because of their relative mechan- 
ical simplicity, high area-to-mass ratio, and compact packaging characteristics. Essential to the design 
and fabrication of inflatable reflectors is the ability to predict the reflector shape upon inflation. When 
neither a deep reflector nor extreme surface precision is required, an attractively simple configuration is 
an initially flat and unstressed circular elastic membrane attached at its perimeter to a stiff ring and sub- 
jected to differential pressure. While uniform loading applied to a constant-thickness membrane will not 
produce the exact paraboloidal reflector shape that is desired in numerous applications, the shape differ- 
ence may be small enough for many purposes. 

This paper, because of inconsistent citations in the literature (e.g., ref. 1) and questions about con- 
vergence, reexamines Hencky’s original analysis (ref. 2) of the large deflections of a clamped, circular 
membrane under uniform pressure. Hencky’s power-series approach is again employed, an algebraic 
error is corrected, more terms are retained to assess convergence, and results are generated for two addi- 
tional values of Poisson’s ratio. 

Also, because Hencky’s problem actually involves uniform lateral loading (i.e., the radial compo- 
nent of pressure on the deformed membrane is neglected), the boundary value problem for true uniform- 
pressure loading is formulated and solved. Results for lateral deflections and membrane stresses from 
both Hencky’s solution (corrected) and the uniform-pressure solution are presented in tabular and 
graphical form, and comparisons are made between the two solutions. 

In addition, although neither problem solution yields an exactly paraboloidal deflection shape, 
nonuniform membrane-thickness distributions that yield such shapes can be found for both loading 
conditions. These distributions are derived in the appendix. 

Symbols 

a radius of membrane 


a 2n> ^2 n* n 2m » w 2 m coefficients in power series 
E modulus of elasticity 

h thickness of membrane 



N dimensionless meridional stress resultant ( N r /(Eh ) in Hencky’s problem, N r /pa in 

uniform-pressure problem) 

A meridional and circumferential stress resultants, respectively 

p uniform loading (lateral in Hencky’s problem, normal in uniform-pressure 

problem) 

q dimensionless loading parameter, palEh 

r, G plane polar coordinates 

u radial displacement 

W dimensionless lateral deflection, w/a 

w lateral deflection 

a constant, (3 + jll)/2 

e r ,£e meridional and circumferential extensional strains, respectively 

(i Poisson’s ratio 

p dimensionless radial coordinate, rla 

Analysis 

The problem of the large deflections of a uniform-thickness, circular isotropic elastic membrane, 
clamped at its boundary without pre- tension and subjected to uniform lateral loading (Hencky’s prob- 
lem), is addressed first. Next, the same configuration under uniform pressure is analyzed. 

Uniform Lateral Loading (Hencky’s Problem) 

The governing equations are 
radial equilibrium 


"e = 



lateral equilibrium 


N dvv = pr 
r dr 2 


( 1 ) 


(2) 


where N r and Nq are, respectively, meridional and circumferential stress resultants, r is the radial coor- 
dinate, w is the lateral deflection, and p is uniform lateral loading; 

stress-strain 


N Q-V- N r = Eht e 
N r -\iN Q = Ehz r 


strain-displacement 



(3) 

(4) 

(5) 


_ du 1 (dw^ 

£r dr + 2\dr J 


( 6 ) 
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where u is the radial displacement and p. is Poisson’s ratio. Equation (2) is the result of one integration 
of the original lateral equilibrium equation (ref. 2) and use of the symmetry condition dw/dr( 0) = 0 
along with the regularity of N r atr = 0. The boundary conditions at the clamped edge are 

w(a) = 0 (7) 


u(a) = 0 (8) 

Combining equations (1) and (3) through (6), and defining dimensionless quantities W = w/a, 
N = N r /(Eh), p = r/a, and q = pa/(Eh), the resulting equations are 


d 

% 




(9) 


, ,dW 1 

N dp m ~2 qp 


( 10 ) 


Substitution of equation (10) into equation (9) gives 




(ii) 


After Hencky’s approach, the following forms for N( p) and VT(p) in equations (10) and (11) are 
assumed: 


1 00 

N(P) = < 12a > 

0 

W( p) = « 1/3 ia 2 „(l-p 2 " +2 ) 02b) 

0 

By these choices, equations (10) and (11) become independent of q 9 so that the power series in 
equations (12a) and (12b), once determined for a specified value of |l, are valid for all q. Also, substitut- 
ing equation (12a) into the radial equilibrium equation gives, for the dimensionless circumferential 
stress resultant, 


Eh 


1 2/3 
4 q 


£(2n+l)b 2n p 2n 


0 


Note that the form assumed for W( p) ensures term-by-term satisfaction of the boundary condition on w y 
equation (7). Then equation (11) becomes 


, 2 , 4,6 

i + &2 P +fe 4P +b 6? + 


-ii 


4(2)fc 2 p + 6(4)& 4 p 3 + 8(6)6 6 p 5 + 10(8)/7 g p 7 + ... = -8p (13) 


3 



Expanding the left side of equation (13) and equating coefficients of like powers of p yields the follow- 
ing relations between b 0 , b 2 , b 4 , b 6 , ..., which can be solved successively for b 2 , b 4 , b 6 , ... in terms of 

bo- 


b 0 b 2 “ 


( 14 ) 


3 b]b A + 2b 0 b j = 0 


( 15 ) 


6blb 6 + 6b Q b 2 b 4 + b 2 (bl + 2b Q b 4 ) = 0 


(16) 


\0b 2 Q bz + \2b Q b 2 b 6 + 3b 4 (b 2 2 + 2b 0 b 4 ) + b 2 (2b Q b 6 + 2b 2 b 4 ) = 0 


(17) 


IS^o^lo + 20 b 0 b 2 b % + bb^bl + 2 b Q b 4 ) + 3b 4 (2b Q b 6 + 2 b 2 b 4 ) 
+ b 2 (b 4 + 2b 0 b^ + 2b 2 b^) = 0 


(18) 


21^ 12 + 30 b Q b 2 b w + 10/?g(Z>2 + 2b Q b 4 ) + 6b 6^ 2b 0 b f, + 2b 2 b 4) 

+ 3b 4 (b 2 + 2b Q b% + 2b 2 b 6 ) + b 2 (2bQb l0 + 2b 2 bg + 2b 4 b 6 ) = 0 


( 19 ) 


In all, the first 10 of these equations have been derived and solved. The results are 

b 2 = ~ 2 

b 0 


( 20 ) 


b 4 = ~—5 
* b 0 


( 21 ) 


b 6 = 


13 

18 (>r 


( 22 ) 


*8 = 


17 


18 b, 


11 


(23) 


^10 “ 


37 


27* 


14 


(24) 


4 



( 25 ) 



] In reference 2, Hencky terminated this sequence with b ^ for which he reported the erroneous value of -407/1 89/? 0 . 
Hencky’s error was corrected by J. D. Campbell in reference 4, where Hencky’s problem also was generalized to include arbi- 
trary initial tension. 

™Mathematica is a registered trademark of Wolfram Research, Inc. 
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For each specified value of p, solution of a suitably truncated version of equation (3 1 ) yields the 
single value of b 0 for which b 0 > 1. (The infinite series in eq. (31) diverges for -1 < b 0 < 1, and 6 0 ^ -1 
implies compressive stresses.) To investigate the convergence of with the number of terms retained, 
a sequence of truncated versions of equation (31) containing from 2 to 1 1 terms can be solved to provide 
a sequence of values of £> 0 . This process has been done for each of the three values of (I considered 
herein. The results with 1 1 terms retained, which differ very little from those with only 6 terms retained 
(as was done in ref. 2), are presented in the following table. 



b 0 

0.2 

1.6827 

0.3 

1.7244 3 

0.4 

1.7769 


a Largely because of an error in the eval- 
uation of 6 I2 , Hencky obtained b 0 = 1.713. 
By apparently truncating equation (12a) 
of the present paper after the b\ 2 term, 
Campbell obtained a value for b 0 of 1 .724, 
which is very close to the value of 1.7244 
found herein. 


With the value of b$ in hand, the truncated power series for both stress resultants and the radial dis- 
placement are completely determined. To obtain the coefficients in the series for W(p), expressions 
(12a) and (12b) are inserted into equation (10), which becomes 


\*> 0 +b 2 p 2 + b 4 p 4 + b 6 p 6 + ... lL 0 + 2a 2p 2 + 3fl 4p 4 + 4 ^6p 6+ f = 1 (32) 


Equating coefficients of like powers of p yields the following relations between the two sets of 
coefficients: 


2b^a 2 &2 a 0 — 0 

3/?qA 4 + 2b 2 a 2 + b<x a Q - 0 

4b 0 a 6 + 3b 2 a A + 2b 4 a 2 + b b a 0 = 0 
5 b 0 a 8 + 4b 2 a 6 + ~ >b 4 a 4 + 2b 6 a 2 + b % a 0 = ® 


(33) 
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Successive solution of equations (33) for a 0 , °2> °4> ••• combined with equations (20) through (29) 
gives 


a 0 ~ u 

1 

a 2 ~ 4 

2b n 


a , i = 


a, = 


9^o 

55 

12b, 


10 


do — 


6b, 


a io - 


13 

0 

205 


108 b 


16 


a 12 - 


'0 

17051 


5292b, 


19 


a 14 ~ 


'0 

2864485 


5080326 


22 

0 


a 16 = 


fl 18 ~ 


103863265 


102876486 
27047 983 


25 

0 


a 20 “ 


1 4696646, 


28 

6 


42367613873 
1 244805 4086q’ 


(34) 


Substitution of these coefficients into equation (12b) gives the dimensionless lateral displacement W(p). 
At this point, truncated power series can be written for all stress resultants and displacements. Tables 
and graphs of stress resultants and displacements are deferred until comparisons can be made with cor- 
responding results from the solution for uniform pressure. 


Uniform Pressure 

All of the equations governing the uniform-pressure problem are identical to those of Hencky’s 
problem (i.e., eqs. (1) through (8)) except for equation (1), the radial equilibrium equation. The new 
radial equilibrium equation can be obtained by application of the principle of virtual work or by 
summing forces on an infinitesimal element of the deformed membrane. By noting that wherever the 
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membrane has nonzero slope the normal pressure has a radial component (neglected in Hencky’s prob- 
lem), the radial equilibrium equation is seen to be 


AT ^ / AT \ d W 

N » = Tp"J-e'T r 

Combining equations (35) and (3) through (6) yields 


n: 


( id N 


dN. 


dr 




p 2 2 dN + ^ 2 2 Ehp 2 r 2 

2 TtH' 2 F rN ' + 8 = ° 


which must be satisfied along with 


N dw = _pr 
r dr 2 


The boundary condition on radial displacement takes the slightly more involved form 


(35) 


(36) 


(37) 


= 0 (38) 

=a 

The nondimensionalization employed in solving Hencky’s problem does not lead here to versions 
of equations (36) and (37) that are independent of q . Nevertheless, the nondimensionalization was used 
initially in the solution process, but was abandoned when inaccuracies were encountered in the numeri- 
cal evaluation of series coefficients. With the definitions p = rta, W = w/a, N = N r /pa,md q = ( pa)/(Eh) y 
equations (36) and (37) become 


[dr 


(rN r )-\lN r - pr 


dW 

dr 


] 


N 


( 2 \ 
1 2 d N _ dN 

o — + 3p 5? 


V dp 
where a = (3 + (J.)/2, and 

Substituting the assumed form 


1 3 dN 2 .. 1 p A 

-p — + ap/V + -— = 0 
2 k dp y 8 q 


dW 1 
N— = --p 


dp 


N (P) = X "2m P 

0 


2m 


(9>0) (39) 


(40) 


(41) 
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into equation (39) and equating coefficients of like powers of p gives the following sequence of rela- 
tions between n 2> n 4> n 6> • • •> which can be solved for n 2 , n 4 , n 6 , ... in terms of kq. 


8no« 2 + an 0 + ^- = 0 
2 2 

24 «o n 4 + ^ n 0 n 2 ( a _ 1 ) n 2 ~ ® 

2 2 
48/10^6 + 48nQ^2 w 4 + 8(2 /Zq« 2 + n 2) n 2 + ( a “ 2)« 4 = 0 

2 2 

80n 0 «g + 96n 0 «2«6 + 24(2« 0 « 4 + n 2 )n 4 + 16 (« 0 n 6 + n 2 n^)n 2 + (ot - 3)n 6 = 0 

2 2 
120rtQrt 10 + 160«Q«2 n 8 + 48(2/io^4 + « 2 ) w 6 + 48(nQ^6 + n 2 n 4^ n 4 


(42) 


+ 8(2«q« 8 + 2n 2 n 6 + n^)n 2 + (a - 4)«g = 0 

2 2 
168n Q n 12 + 240n 0 n 2 Ai 10 + 80(2« 0 n 4 + n 2 )n g + 96(n 0 n 6 + « 2 n 4 )« 6 


+ 24(2n 0 n 8 + 2« 2 n 6 + n 4 )n 4 + 16(« 0 n ]0 + n 2 n g + n 4 n 6 )n 2 + (a-5)n 10 = 0 


The first eight of these equations have been derived and solved, although only the first six are shown in 
equations (42) because of rapidly increasing algebraic complexity. For the same reason, only the first 
four solutions are shown here. 


n 2 


n 


4 


n 6 

n 8 


(1 + 8a«Q<7) 

64hq<? 

(1 + 8 an 0 < 7 )( 1 + 4 n Q q + 4a n Q q) 

6144«q^ 2 

2 2 2 2 2 2 
(1 + 8a« 0 ^)(13 + 96n 0 q + \2San Q q + 128n 0 ^ + 576a« 0 + 256a n Q q ) 

4718592n^ 3 

-(1 + 8a« 0 ^)(39 + 366n G q + 670a n Q q + 288/ig# + 960 n^q 2 + 460&an^q 2 
+ 3584a 2 «g<7 2 + 4608nJ$ 2 + 768n 3 $ 3 + 6272a«g? 3 + 14208a 2 n^ 3 
+ 5632a 3 n 3 0 q 3 + 36864a n 5 Q q 3 - 1 8432a 2 ng<7 3 )/754 974720nJ V 


(43) 


In the development of equations (42) and (43), the Mathematica system was again used extensively. 
Substitution of equations (43) into equation (41) yields the power series for N(p) in terms of the yet- 
undetermined coefficient /ig that, in this problem, depends on both p and q. However, before the 
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remaining boundary condition, equation (38), can be applied to evaluate n 0 , equation (40) must be used 
to express the coefficients in the assumed series 


W(P) = £ vv 2n (l-p 2 " +2 ) (44) 

0 

in terms of n 0 . Substituting equations (41) and (44) into equation (40) and equating coefficients of like 
powers of p gives the system of simultaneous equations 


= \ 

2 n 0 w 2 + n 2 w 0 = ® 

3«qW 4 + 2n 2 w 2 + n 4 w o = ® 

4 "0^6 + 3n 2 w 4 + 2n 4 w 2 + n 6 w 0 ~ 0 
SriQWft + ^ n 2 w e + ^ n 4 w 4 + ^ n 6 w 2 + ^8^0 = ^ 


(45) 


from which vv 0 , vv 2 , vv 4 , ... can be found in terms of n$. For the sake of brevity, only the first five solu- 
tions are shown: 

w 0 = l/4n 0 

4 

w 2 = (1 + &anQq)/5\2n 0 q 

1 2 

w 4 = ( 1 + 8 cx/1q#)(5 + 8/1q^ + 32attQ<7)/147456«Q<7 
w 6 = ( 1 + 8aM Q ^)(55 + 192n 0 ^ + 704an 0 ^ + 128 «q^ 2 + 1344a/iQ<j 2 
+ 2176a 2 «o<7 2 )/75497472>4V 

w 8 = (1 +8an 0 ^)(259+ 1366n 0 <? + 5030coi 0 ? + 288/1^ + \920n\q 1 (J6) 

+ 19008an^ 2 + 31 lAAanlq 1 + 4608 n Q q 2 + 768 h^q 3 
+ 13952an^ 3 + 65408a 2 n^ 3 + 64512a 3 /i 3 ^ 3 + 36864an^ 3 
- 18432a 2 no< ? 3 )/15099494400/i ( ! | V 


Now n 0 can be evaluated through satisfaction of the boundary condition on radial displacement, 
equation (38), which has the dimensionless form 



( 47 ) 
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Substitution of equations (41) and (44) into equation (47) yields 

(^} i 0 + (^) , 2 + (^} 1 4 + (^} i 6 + --- + vv 0 +2m, 2 + 3w, 4 + 4 vv 6 + --- =0 < 48 > 

in which w 0 , vv 2 , w 4 , ... are given in terms of n 0 by equations (46). 

Substitution of equations (43), giving n 2j n 6 , ... in terms of tt 0 , and equations (46), giving w 0 , 

w 2 , w 4 , ... in terms of mq, into equation (48) yields an equation for the lead coefficient After specifi- 
cation of values for |i and q y various truncations of this equation can be solved until a satisfactorily con- 
verged value of hq is obtained. This value of can then be used along with equations (43) to write the 
explicit truncated series for N(p). Similarly, equations (46) can be inserted into equation (44) to produce 
the series for W(p). With these two series in hand, expressions for all other variables of interest can be 
generated. 

Results and Discussion 

Lateral Deflections 

The variables of greatest interest in the design of membrane reflectors are the lateral deflection w 
and the meridional stress resultant N r . In figure 1, the dimensionless lateral deflection W = wla for 
Hencky’s problem is plotted as a function of p = rla , the dimensionless radial coordinate, for q = 0.001 , 
0.01, and 0.1 and 1 1 = 0.2, 0.3, and 0.4. 

Note that, in Hencky’s problem, the family of curves for a fixed value of q can be used to produce 
curves for any other value of q by appropriate adjustment of the vertical scale. The same is not true, 
however, in the solution to the uniform-pressure problem. Figure 2 contains plots of the lateral deflec- 
tion from the uniform-pressure problem for the same set of values of p and q shown in figure 1 . In both 
solutions, as deflections increase with q , so do the Poisson’ s-ratio-related differences among them, with 
larger values of p corresponding to smaller deflections. 

To illustrate the effect on lateral deflection of retaining the radial component of the pressure, results 
from figures 1 and 2, for p = 0.3 only, are plotted together in figure 3 for q- 0.001, 0.01, and 0.1. For 
q - 0.001, the deflections for the two loading cases are nearly identical, as would be expected, because 
the radial component of the pressure is proportional to dw/dr , which is everywhere small when q is 
small. In fact, even for q = 0.01, the deflection differences are still very small although they could have 
some significance for application to high-precision surfaces for electromagnetic reflectors. For q - 0.1, 
a value that corresponds to very high membrane stresses, a basic difference between the deflection 
shapes is quite apparent. The uniform-pressure loading causes a more nearly spherical shape, one that is 
even farther removed from the ideal paraboloid than the shape predicted by the solution to Hencky’s 
problem. Thus, the present results indicate that a reflector design based on the solution to Hencky’s 
problem would likely overestimate electromagnetic performance. With regard to achieving the ideal 
deflection shape, nonuniform membrane-thickness distributions that lead to paraboloidal middle- 
surface deflection shapes can be found in both problems. Details of the analyses are presented in the 
appendix. 

Another useful way of viewing the results is by envisioning the loaded membrane as a nonlinear 
spring. Figure 4 contains plots of the dimensionless loading parameter q as a function of the dimension- 
less center deflection w(0)/a from each solution for the three values of Poisson’s ratio. Both solutions 
depict systems that are stiffening with increasing load, as evidenced by the increasing slopes, especially 
at the larger values of q , with the normal pressure solution representing the stiffer system. 
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Stress Resultants 


Figure 5 contains plots of the dimensionless meridional stress resultant N r /(Eh) for q- 0.001, 0.01, 
and 0.1, with p = 0.2, 0.3, and 0.4. For q = 0.001, the curves for Hencky’s problem and for the normal- 
pressure case are indistinguishable on the figure. For q = 0.01, differences between the two solutions 
still appear to be quite small, particularly at the scale used for the figure. The differences become more 
pronounced as q continues to increase until, for q = 0.1, the uniform-pressure stresses substantially 
exceed those of Hencky’s problem at the center, while the opposite is true at the boundary. The follow- 
ing table contains a list of values of the maximum dimensionless stress resultant N r (0)/(Eh) from both 
solutions as a function of q for all three values of p. 


Q 

Uniform pressure 

Hencky 

p = 0.2 

p = 0.3 

*P 

II 

O 

p = 0.2 

T= 

li 

O 

u> 

p = 0.4 

0.001 

0.00425 

0.00435 

0.00448 

0.00421 

0.00431 

0.00444 

.002 

.00677 

.00693 

.00714 

.00668 

.00684 

.00705 

.003 

.00890 

.00911 

■ 


.00897 

.00924 

.005 

.0126 

.0129 

■ 


.0126 

.0130 

.008 

.0173 

.0177 

.0183 

.0168 

.0172 

.0178 

.010 

.0202 


.0213 

.0195 

.0200 

.0206 

.020 

.0326 


.0343 

.0310 

.0318 

.0327 

.030 

.0433 

.0443 

.0456 

.0406 

.0416 

.0429 

.050 

.0622 

.0637 

.0655 

.0571 

.0585 

.0603 

.100 

.1032 

.1057 

.1086 

.0906 

.0929 

.0957 


Thus, while it appears that for all q > 0 the solution to Hencky’s problem understates the maximum 
membrane stress resultant (;V r (0), which is equal to Nq(0) in either solution), the discrepancy is substan- 
tial only for large q , in which case the stress resultants predicted by either theory are very large. 

Concluding Remarks 

The stresses and large deflections of a clamped, circular membrane have been analyzed for two 
types of loading: uniform lateral loading (Hencky’s problem) and uniform-pressure loading. Results for 
the two types of loading have been shown to agree quite closely when the membrane is lightly loaded 
and to diverge slowly as load intensity is increased. The deflection shapes from the uniform-pressure 
solution are more nearly spherical than those from the solution to Hencky’s problem, especially at the 
higher load intensities. This result suggests that a near-paraboloidal inflatable reflector designed on the 
basis of the more convenient Hencky solution would overestimate electromagnetic performance 
because of shape-accuracy shortcomings. 

Both problems were found to yield closed-form nonuniform membrane-thickness distributions that 
produce exact paraboloidal middle-surface deflection shapes. The fabrication of such variable-thickness 
membranes appears to pose a daunting challenge. 


NASA Langley Research Center 
Hampton, VA 23681-2199 
March 11, 1997 
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Appendix 

Membrane-Thickness Distributions That Yield Parabolic Deflection Shapes 

To find thickness distributions that yield parabolic middle-surface deflections under uniform load- 
ing, the governing equations are approached under the assumptions that w(r) is a prescribed parabola 
and the thickness h(r ) is to be determined. The variable-thickness version of Hencky’s problem is 
addressed first, followed by the uniform-pressure problem. 


Hencky’s Problem 

The governing equations are repeated here for convenient reference. 


= J r W. r> 

(Al) 

i N i - -EL 

r dr 2 

(A2) 

Ehz r = N r -\iNQ 

(A3) 

Eh H = N B -VLN r 

(A4) 

_ du \(dw Y 

£ r dr + 2\dr J 

(A5) 

u 

£ e = ; 

(A6) 


A parabolic deflection shape is assumed. 


w(r) 


f 

w 0 + W 1 

V 


2\ 

r 

~2 
a ) 


(A7) 


where w 0 is an unspecified rigid-body translation and W is the strain-related maximum deflection, so 
that 


dw 

~dr 


2 Wr 


2 

a 


Then equations (Al) and (A2) yield 


a/ - A r - 

N r~ N Q~ 4W 


Combining this result with equations (A3) through (A6) gives 


Eh 


du W 2 r 2 ^ 

* a 4 , 


= (1 -P-) 


pa_ 

4 W 


and 


= (l-(t)p fl 2 r 

4 WEh 


(A8) 


(A9) 


(A 10) 


(All) 
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Substitution of equation (All) into (A 10) and some simplification yield 


i(l] = 8EW 3 r 

dr ^> = ~0- i l)pa 6 

After integration, with the definition h( 0) = /iq, the solution is found to be 


r 4Eh n W r \ 

1 2 6 

- (1 - 1 *)pa - 

which, on physical grounds, requires satisfaction of the inequality 


(A 12) 


(A 13) 


4Eh 0 W 3 

(1 -p)pu 4 


(A 14) 


Note that, for sufficiently small values of the quotient in the denominator of equation (A 14), the 
radial thickness variation is very nearly parabolic. The focal length F of the parabolic deflection shape is 
given by F = a 2 /4W, so that W = a/8(F/D ), where D = 2a is the diameter of the circular reflector. Then, 
in terms of dimensionless ratios, an alternate form for equation (A13) is 


h(r) 


1 E h _Q {rlaf^ 
128(1 - p) p a 


Also, equation (A9) can be rewritten as 


N r = = P F = Vo 


(A15) 


(A 16) 


in which Co is the (maximum) membrane stress at r = 0, where h(r ) is minimum. Thus, with the focal 
length F prescribed and an acceptable working value chosen for a 0 , equation (A16) enables a family of 
choices for p and h$ that, along with values for D and material properties E and p, can be inserted into 
equation (A15) to produce a family of candidate thickness profiles from which to choose. Clearly, the 
larger the value of FID (i.e., the shallower the reflector) the greater the freedom in assigning values to 
the remaining parameters. 

Note that this solution cannot be made consistent with the original boundary condition u(a) = 0. 
From equations (A4), (A6), and (A9), the radial displacement takes the form 


u(r) 


A 2 
1 - p pa r 

E 4Wh(r) 


(A 17) 


from which it is clear that u(r) must be positive for all r > 0. Thus, in the variable-thickness version of 
Hencky’s problem, achieving the parabolic deflection shape would also require an appropriate uniform 
expansion of the circular boundary. 


Uniform-Pressure Problem 

The governing equations for this problem are equations (A2) through (A6) along with 

£<r/V,)-N e = pr^ (A 18) 
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in place of equation (Al). Assumption of the same parabolic displacement function, equation (A7), 
leads to 


N =££_ 
r 4w 


(A19) 


and 


N e = m +2pwr - 


(A20) 


Insertion of equations (A 19) and (A20) into equations (A3) and (A4) followed by substitution of the 
results into equations (A5) and (A6) leads to 


Eh 


r du 

dr a A J 


V 


2 2 
a 


(A21) 


Eh - r = ^-^ + 2pW -l 

a 


(A22) 


Eliminating u(r ) between equations (A21) and (A22) gives 


4 

dr 1 h 


2 u/ 3_ i 

1 4W 2 


2 2~l 

(l_M)g£L_2 ^ Wr 
U P) 4W 2 


a 


-2 


2 2 

EW r 


(A23) 


Differentiation followed by simplification, where possible, leads to a linear first-order differential equa- 
tion in 1 !h(r) with variable coefficients. Its solution with the help of an integrating factor yields 


h{r) = A 0 - 


1 + 


sirV ' l(3+tA)/2 


(i-n)« J 


i - 


EWh 


0 


(3 + p )pa 


1 + 


2 2 

8vr> 

(1 -|x)a 4 


(3+p)/2 


(A24) 


- 1 


where, again, h 0 = h(0). Clearly, h(r) has its minimum value at r = 0. 

An alternate form of equation (A24) in terms of the dimensionless ratios employed earlier is 


h(r ) = h Q - 


1 + 


(r/aY 


-,(3+p)/2 


8(1 -\i)(F/dY 


1 E h o 1 


8(3 + ji) p a ( F/ £)) 2 1[ 8(1 -p)(F/D) 2 J 


1 + 


(r/aY 


n(3+g)/2 


(A25) 


- 1 


This equation coupled with 


N r (0) = W e < 0) = PL 
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or, equivalently. 


pF = a 0 h Q (A26) 

can be used along with combinations of material and geometric properties, some of which are likely to 
be mandated by a particular application, to produce various candidate thickness distributions from 
which to choose. As in the preceding problem, achieving the parabolic deflection shape would require 
appropriate expansion of the circular boundary, which is apparent from equation (A22). 

The fabrication of a membrane with a thickness distribution mandated by either of these theories 
may be beyond the current state of film manufacturing technology. Of the two families of thickness 
distributions, the one appropriate to Hencky’s problem may be nearer to achievable, possibly by use of 
a rotating, ultraflat film-fabrication table, because it is very nearly parabolic when the thickness varia- 
tion relative to the minimum thickness is not too great. 
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Figure 2. Lateral deflection due to uniform-pressure loading. 
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Figure 3. Lateral deflections for |i = 0.3; q = 0.001, 0.01, and 0.1. 
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Figure 5. Meridional stress resultant distributions for q = 0.001 , 0.01 , and 0. 1 . 
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